Quality factor of a matter-wave beam. 
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Imperfections in dilute atomic beams propagating in the paraxial regime and in potentials of cylin- 
drical symmetry have been characterized experimentally through the measurement of a parameter 
analogous to a beam quality factor [Riou et al., Phys. Rev. Lett. 96, 070404 (2006)]. We propose 
a generalization of this parameter, which is suitable to describe dilute matter waves propagating 
beyond the paraxial regime and in fully general linear atom-optical systems. The presented quality 
factor shows that the atomic beam symmetry can be traded for a better transverse collimation. 

PACS numbers: 03.75.Pp, 41.85.Ew, 31.15.xh 



I. INTRODUCTION 



There is a strong analogy between the propagation of 
light and matter waves 1:]. It is manifest in the paraxial 
wave equation for the mode U of an electromagnetic field 
E(x,y,z) — U(x,y,z)e % ( kz ~ ut 'e which takes the form of 
a bidimcnsional Schrodinger equation 



2ik 



dU(x,y,z) 
dz 



dx 2 dy 2 



U(x, y, z)+F(x, y)U{x, y, z) 
(1) 

with a function F(x,y) depending on the local refraction 
index. A spectacular achievement proving the similarity 
of the atomic and light fields was the realization of a 
quasicontinuous matter beam analog to a laser It 
has been shown recently that the extraction of a quasi- 
continuous atom laser from a Bose-Einstein condensate 
by a weak RF [lB - [l3l ] or Raman [3] outcoupling involves 
the propagation through a discontinuous potential which 
degrades the beam collimation. To quantify this effect, 
Riou et al. fTo| introduced a parameter analogous to 
a quality factor. It allows one to estimate the the 
transverse width of a cylindrical atomic beam falling in 
a uniform gravitational field and in the paraxial regime, 
i.e., with atoms strongly accelerated along one privileged 
direction. 

However, an important difference between light 
and matter waves is that the latter generally do not 
propagate in this paraxial regime. This is typically 
the case for a condensate expanding after a sudden 
trap shut down, which can be viewed as a pulsed atom 
laser. Even for a weakly outcoupled atom laser, the 
paraxial approximation is often valid only after a certain 
propagation time [lOj during which the beam must be 
treated in a more general framework. A comparison 
between the usual Schrodinger equation and the paraxial 
wave equation (JTJ shows that for these matter waves, 
time plays the role devoted to the axis of propagation 
Oz in paraxial optics, and that the transverse space is 
three dimensional and spanned by {e x ,e y ,e z } instead 
of {e x ,e y }. 

To extend the concept of matter beam quality beyond 
the paraxial regime, it is necessary to reconsider its 



physical significance in optics. For cylindrical light 
beams, a quality factor compares the divergence of a 
non-ideal beam to a standard set by a perfect Gaussian 
mode. It conveys the idea of degradation in the collima- 
tion or, equivalently, that of diminution in phase-space 
density. A relevant quality parameter should thus be 
left invariant during the propagation in perfect linear 
optical systems which preserve the beam collimation. 

Such systems are modelled by a finite set of ABCD 
matrices which describe linear input-output relations in 
the light beam phase-space [l5j]. Their atomic coun- 
terpart is the quantum evolution under an Hamiltonian 
quadratic in position and impulsion, which involves simi- 
lar relations with 3x3 ABCD matrices [H, [T3|- The 
phase space evolution of the atomic beam then also 
amounts to a time-dependent map in the arguments of 
the Wigner distribution [l[ 

W (r, p,t) = W[ D(r - f ) - — B{p - m4>) , 



-mC(r - £) + A(p - m<j>) , t 



The ~ stands for the transposition, the vectors r,p are 
the position and impulsion, the matrices A, B, C, D 
and vectors £,$6 are time-dependent parameters deter- 
mined in Ref.[16]. Thanks to the unimodularity of the 
ABCD matrices, the propagation under those quadratic 
Hamiltonians preserves the phase-space density and as 
such do not alter the atomic beam collimation. In this 
view, these Hamiltonians can be considered as perfect 
aberrationless atom-optical systems and should thus 
leave invariant a well-defined quality factor. 

This is not the case for the parameter used so far to 
characterize the quality of matter beams [l(l 13 • This 
is disturbing, since its current definition would some- 
times lead to ignoring possible improvements of the atom 
laser collimation through non cylindrical atom-optical el- 
ements such as astigmatic atomic lenses 18]. It seems 
thus useful to propose a new quality factor for atom lasers 
which respects this invariance requirement. This is the 
main purpose of this paper. 
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II. QUALITY FACTOR OF A CYLINDRICAL 
BEAM 



Let us first reconsider the definition of a quality factor 
for a cylindrical light beam, before extending this con- 
cept to atom optics. Partially coherent light beams can 
be described by means of a first-order field correlation 
function T relating the field amplitude at different points 
of planes transverse to the direction of propagation O z . 
For a cylindrical beam, only a single transverse direction 
O x needs to be considered in the correlation function 

T(x 1 ,x 2 , z) = (E(xi, z) E*(x 2 , z)) . 

The Wigner transform of this function provides a phase 
space picture of the beam 



W(x,k x ,z) 



J dx >r(x + ^,x-^,z)e-* k * x ' , 



which can be used to define moments in the transverse 
position and wave vectors 

(m(x, k x ,z)) = J Ax dk x m(x, k x ,z) W(x, k x ,z) . 

Because the beam is cylindrical and transverse 
(a;) = (k x ) = 0. The moments Ax = yj (x 2 ) and 
Afc-r = \J (k 2 ) are, respectively, the transverse squared 
width and wave- vector dispersions of the beam. 

Optimal collimation is achieved with a perfectly coher- 
ent Gaussian mode, for which the position and direction 
moments verify at the waist 



^ X L, Gaussian ^^L,Gaussi 



1 



Partially coherent cylindrical light beams have a lessen 
collimation quantitatively estimated thanks to a quality 
factor M 2 This parameter can be expressed as a 

combination of moments left invariant during the propa- 
gation in aberrationless cylindrical optical systems 



M 2 



(2) 



At the beam waist, this expression reduces to the product 
of transverse squared width and divergence which has a 
clear interpretation in terms of phase-space dispersion 



M2 



A a;,,, Afc^. 



(3) 



This is indeed the definition which has been adopted to 
express the quality of a matter beam [l(| E3] > once taken 
into account De Broglie relation between momentum and 
wave vector p = Kk 



: Aa;L Ap x 



(4) 



This definition is consistent if one considers only 
matter-waves and atom-optical systems with cylindrical 
symmetry, and a propagation in the paraxial regime: 
the parameter M 2 is then unaltered by the propagation 
in aberrationless systems. 

However, these conditions are far too restrictive to 
characterize the matter beams in most experiments, 
which generally involve potentials without any symme- 
try and a matter-wave propagation outside the paraxial 
regime. A three-dimensional quality factor M3, invariant 
under the whole set of possible ABCD matrices, is then 
required. Following previous treatments of optical and 
quantum mechanical invariants [20I422T ] , we now examine 
how to properly define this parameter. 



III. QUADRATIC PROPAGATION INVARIANT 
FOR MATTER BEAMS 

Our goal is to replace the parameter Mi by an ABCD- 
invariant combination of moments which still provides 
an insight on the matter beam divergence. In the unidi- 
mensional case, the price to pay to express the product 
Ax| w Ak x \ w ^ as an invariant expression © was to 
introduce an additional moment (x k x ) mixed in wave 
vector and position. This term can be viewed as an el- 
ement of a variance matrix. Such matrices have been 
introduced in optics [20|, and they can also characterize 
a partially coherent atomic beam [23[ : 



V = 



A rr A r 



The matrices A a b are 



A 



ab 



(a x b x ) (a x b y ) (a x b z ) 
{a y b x ) {ayb y ) (a v b z ) 
(a z b x ) (a z b y ) (a z b z ) 



with the vectors a, b = r, p/m. We now derive the vari- 
ance evolution when the atomic field propagates under an 
Hamiltonian which is quadratic in position and impulsion 



H(r,p) 



P0{t)P 
2m 



ra(t)p ?7(i)r— mg(t)-f +f (t)-p . 

(5) 



a(t), f3(t) and j(t) are 3x3 matrices; f(t) and g(t) are 
three dimensional vectors. It follows from the equations 
of motion for the position and impulsion operators in the 
Heisenberg picture that the variance matrix V satisfies 



dt 



VT 



r(t) 



a(t) m \ 
lit) a(t) J 



The integration of this last equation involves the time- 
ordering operator T and the atom-optical ABCD ma- 



3 



trix 124 



M(t,t ) 



( A(t,t ) B(t,t ) 
\C(t,t ) D(t,t Q ) 



Its expression in terms of position and impulsion mo- 
menta is 

M s = ^2 [( x2 )(pl) - ( x Px) 2 + (y 2 )(p 2 y) - (ypy) 2 



= T exp 



- / dt' 

to 



a(t') fi(t<) 
lit') a(t>) 



3ft 2 

+ (z 2 )(pt) - (zpz) 2 + 2({xy){p x p y ) - (xp v ) (yp x )) 



+ 2({xz){p x p z ) - {xp z ){zp x )) + 2{(yz){p y p z 



The evolution of the variance matrix in the quadratic 
Hamiltonian ([5]) is then similar to the transformation 
laws in optics 



(yPz)(zpy))] ■ 
(10) 



V(t) = M{t,t )V{t )M{t,t ) 



(6) 



This expression can be rewritten in a more compact form 
as a sum of moments along single and multiple directions: 



The most straightforward attempt to generalize the uni- 
dimensional quality factor M 2 to three dimensions would 
consist in considering norms and scalar products in 
Eq-© instead of single coordinates 

^1 {x 2 + y 2 + z 2 ) {p 2 x + p 2 y + pi) - (r ■ p) 2 . 

At an instant satisfying (r -p) =0, the multidimensional 
equivalent of a "waist" , this expression would express 
the greater phase space occupied by the matter beam. 
Unfortunately, this quantity is not left invariant in the 
propagation under a general quadratic Hamiltonian. 

Nonetheless, it is possible to define a quality factor 
which respects the invariance requirement and still pro- 
vides a meaningful insight into the beam phase space. 
We use the fact that the ABCD matrices associated with 
matter wave propagation are symplectic, i.e., they verify 
at all times the following relation 



M, = 



^ [Qx ~\~ Qy ~T" Qz A X y A XZ Ayz\ 



Q q 



h 2 



(11) 



The quantities Q v correspond to the fourth power of the 
previous matter- wave quality factor Mi considered along 
the three spatial directions, while the terms A v v > reflect 
correlations between different directions. If one chooses 
the coordinate system along the beam principal axis, the 
position momenta satisfy {i] rj') = for r/ =/= r/' , yielding 
a simpler expression for A v v i : 



A , — 
v v — 



These two families of parameters have different physical 
significance and obey different constraints. The parame- 
ters Q v are a manifestation of the beam divergence and 



M(t, t ) — L M{t, t ) L L — M(t, t ) L M(t, t ) , are bounded below by the Heisenberg principle 



with the matrix 



(7) 



L = i 



-1 

1 



This symplectic structure can indeed be used to generate 
a family of invariants [2fj| - |22l | . The lowest order of this 
invariant family indeed generalizes the definition of the 
current matter beam quality factor 



M£(t) 



Am 



Tr [V(t) L V(t) L] 



(8) 



The constant is adjusted to yield Ma — 1 for a Gaussian 
matter wave. Proving the invariance of this parameter is 
simple and identical to optics [20| . Combining relation 
© describing the variance matrix propagation with the 
cyclic property of the trace, the definition (|8]) can be 
recast as 

M 3 4 (t) = |^Tr [v(to)M(t,to)LM(t,to)V(t )M(t,t )LM(t,io)] 

The symplectic relation ([7]) then gives the desired invari- 
ance of Ms{t): 



M 3 4 (i) 



Am 

'm 2 



Tr[V(to)LV{t )L]=Ml{t ). (9) 



, >1. 

The terms A v ^ reflect the beam asymmetry and can 
be of either sign, they cancel for spherical clouds. As 
previously announced, a Gaussian matter wave satisfies 
M 3 = 1. 

This atomic beam quality factor thus reflects the de- 
parture from a fundamental Gaussian mode: a beam with 
a quality factor M3 3> 1 needs to be expanded onto sev- 
eral modes, which is likely to degrade the fringe pattern 
in an atomic interference experiment or to add additional 
noise in the atomic beam [251 ]. In principle, if these modes 
were put into a fully coherent and accurately controlled 
superposition, a high value of the quality factor would 
merely induce an additional complexity in the fringe pat- 
tern. This ideal situation is, however, not encountered in 
practice, since the coefficients of the mode decomposition 
are generally not accessible. For most precision interfer- 
ometric measurement, a simple TEMoo mode is indeed 
highly preferable. In the LIGO experiment, the beam 
quality factor of the laser involved is maintained at a 
value M ~ 1 thanks to a mode-cleaning step which fil- 
ters out high order modes [2^,[23]- This quality factor re- 
quirement, reflecting the control on the beam structure, 
should also apply to accurate interferometers involving 
atom lasers. 
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It seems in general more advantageous to minimize the 
quantities Q n which are directly related to the beam colli- 
mation. In this respect, the fact that the invariant quan- 
tity be M3 instead of the coefficients Q v has a practical 
interest. Considering for instance a cylindrical beam with 
given Q x = Q y = Q r = Mf and Q z , it is possible to re- 
duce Q r while keeping the quality factor M3 constant by 
increasing Q z or by generating negative asymmetrical pa- 
rameters A v n i [28[ . One can thus trade some longitudinal 
collimation or cylindrical symmetry for a better trans- 
verse collimation. The corresponding ABCD transfor- 
mation could be implemented by astigmatic atom-optical 
lenses such as electromagnetic waves with an asymmet- 
ric wave- front 11811 . This possibility of improvement is ig- 
nored in Ref . [lOl fl4| , which assume that the transverse 
divergence measured at the waist (i.e. M±) is an upper 
bound for the atom laser collimation in the subsequent 
propagation. 



of an atom laser. Higher order invariants propor- 
tional to I k = (f ) 2fc Tr [V(t) L V(t) L] h could also be 
considered to describe the atomic beam, but they do 
not provide the same insight into the phase-space density. 

Note: Since the publication of this paper, a different 
nonlinear matter wave quality factor has been pro- 
posed [29] , which is preserved in presence of mean-field 
interactions. However, this nonlinear parameter is 
suitable only for uniform and paraxial atomic beams 
propagating in transverse potentials of cylindrical 
symmetry. The generalizations of the matter wave 
quality factor presented in this paper and in Ref. (29j 
are indeed complementary: they enable one to apply 
the concept of beam quality beyond the paraxial and 
beyond the linear regime of propagation respectively. 
The non-paraxial and the non-linear matter wave quality 
factors can thus be fruitfully applied in different contexts. 



IV. CONCLUSION 



A quality factor has been defined for matter waves 
[M 3 , Eq. fTU]) ]. which addresses the general propagation 
of an atom laser. It generalizes the currently adopted 
beam quality factor [Mi, Eq. flU)], indeed only appro- 
priate to describe the propagation of matter-waves in 
the paraxial regime and in cylindrical potentials, and 
which can lead to underestimate the optimal collimation 
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